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ON  THE  LONG  PERIOD  LUNI-SOLAR 
EFFECT  IN  THE  MOTION  OF  AN 
ARTIFICIAL  SATELLITE 


by 

Peter  Musen 
Goddard  Space  Flight  Center 


SUMMARY 


Two  systems  o£  formulas  are  presented  for  the  determination  of  the 
long  period  perturbations  caused  by  the  Sun  and  the  Moon  in  the  motion  of 
an  artificial  satellite.  The  first  system  can  be  used  to  determine  the  lunar 
effect  for  all  satellites.  The  second  method  is  more  convenient  for  finding 
the  lunar  effect  for  close  satellites  and  the  solar  effect  for  all  satellites. 

Knowledge  of  these  effects  is  essential  for  determining  the  stability  of 
the  satellite  orbit.  The  basic  equations  of  both  systems  are  arranged  in  a 
form  which  permits  the  use  of  numerical  integration.  The  two  theories  are 
more  accurate  and  more  adaptable  to  the  use  of  electronic  machines  than 
the  analytical  developments  obtained  previously. 
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ON  THE  LONG  PERIOD  LUNI-SOLAR 


EFFECT  IN  THE  MOTION  OF  AN 
ARTIFICIAL  SATELLITE 

by 

Peter  Musen 

Goddard  Space  Flight  Center 

INTRODUCTION 


Two  systems  of  formulas  are  presented  herein  for  determining  lunar 
and  solar  long  period  effects  of  the  first  order  in  the  motion  of  an  artificial 
satellite.  The  first  method  is  based  on  the  theory  originally  developed  by 
Gauss  (Reference  1)  for  a  numerical  treatment  of  the  very  long  period  ef¬ 
fects  in  planetary  motion,  which  was  found  to  be  applicable  to  the  case  of 
artificial  satellites.  The  second  method  is  valid  for  close  satellites  and  is 
based  on  the  development  of  the  disturbing  function  in  terms  of  Legendre 
polynomials. 

Knowledge  of  these  long  period  effects  is  essential  in  determining  the 
stability  of  the  orbit  and  the  lifetime  of  the  satellite.  As  an  example,  Kozai 
(Reference  2)  found  that  the  solar  and  lunar  perturbations  have  shortened 
the  lifetime  of  Explorer  VI  (  1959  Delta)  by  a  factor  of  ten.  The  long  period 
lunar  effect  can  be  deduced  by  averaging  the  perturbations  with  respect  to 
the  mean  anomaly  of  the  Moon  and  with  respect  to  the  mean  anomaly  of  the 
satellite.  The  first  averaging  process  is  performed  analytically,  the  second 
process  numerically.  The  long  period  solar  effect  is  obtained  by  averaging 
the  perturbations  with  respect  to  the  mean  anomaly  of  the  satellite  only. 

Both  methods  result  in  the  numerical  integration  of  the  equations  for  the 
variation  of  elements. 
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The  interval  of  integration  depends  upon  the  proximity  of  the  satellite  to 
the  earth  and  upon  the  secular  changes  of  the  node  and  perigee  produced  by 
the  earth's  oblateness.  In  a  normal  case  the  interval  of  integration  will  be 
of  the  order  of  several  days.  For  more  distant  satellites  it  can  be  of  the 
order  of  one  month  or  more.  The  choice  of  such  a  large  interval  is  impossi¬ 
ble  if  Cowell's  method  of  integration  in  the  rectangular  coordinates  is  used. 

For  more  distant  satellites  it  was  found  that  the  development  of  the  dis¬ 
turbing  function  into  series  of  Legendre  polynomials  converges  so  slowly 
that  it  is  impossible  to  include  all  important  long  period  terms  in  an  analyt¬ 
ical  development.  A  large  orbital  inclination  might  also  contribute  substan¬ 
tially  to  the  slowness  of  the  convergence.  In  such  a  case,  as  for  example 
a  =  10  earth  radii,  e  =  0.8,  the  analytical  development  obtained  by  Musen  and 
Bailie  (Reference  3)  becomes  incomplete.  These  circumstances  gave  rise  to 
the  investigation  of  the  possible  use  of  the  Gaussian  method  for  determining 
the  lunar  long  period  effects  in  the  motion  of  an  artificial  satellite. 

Halpnen's  form  (Reference  4,  with  corrections  by  Goriachev,  Reference 
5 )  of  the  Gaussian  methodwas  found  to  be  the  most  convenient ,  partly  because 
it  is  very  adaptable  to  the  use  of  electronic  computers.  Some  necessary 
modifications  were  made  to  include  the  basic  perturbations  of  the  motion  of 
the  Moon.  Also,  the  Goursat  transformation  (Reference  6)  and  the  Euler 
summability  process  were  employed  to  speed  up  the  convergence  of  the 
hypergeometric  series.  This  method  is  valid  for  all  values  of  e,  i  and  a/a'. 

The  development  in  terms  of  Legendre  polynomials  can  be  used  if  the 
lunar  perturbations  for  close  satellites  or  the  solar  perturbations  for  all 
satellites  are  to  be  determined.  The  basic  equations  of  these  two  cases  are 
arranged  in  the  form  which  permits  the  use  of  the  numerical  integration  and 
includes  all  the  significant  long  period  terms.  Theoretically,  these  equations 
are  completely  equivalent  to  the  extensive  analytical  developments  obtained 
by  Kozai  (Reference  2)  and  by  Musen  and  Bailie  (Reference  3)  but  they  have 
a  more  compact  and  symmetrical  form  and  are  more  adaptable  to  the  'use  of 
electronic  machines.  Thus  the  formulas  given  here  represent,  from  a  prac¬ 
tical  point  of  view,  a  substantial  improvement  over  the  previous  methods. 


LUNAR  DISTURBING  FORCE  IN  THE  MOTION  OF  A 
DISTANT  ARTIFICIAL  SATELLITE 

Let  gQ,  CO,  n,  i,  e,  a,  n  be  the  osculating  elements  of  the  artificial 
satellite  referred  to  the  earth's  equator,  and  gp ,  co'  ,  Q'  ,  i',  e'  ,  a'  ,  n'  be 
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the  elements  of  the  Moon  referred  to  the  fixed  ecliptic.  We  assume  that 
there  is  no  sharp  commensurability  between  the  mean  motions  n  and  n'.  It 
is  sufficient  to  take  into  account  the  secular  changes  in  the  oi'  ,  H'  and  to 
neglect  the  periodic  effects  and  the  influence  of  the  precession.  Let  P  be 
the  unit  vector  directed  from  the  center  of  the  earth  to  the  perigee  of  the 
orbit  of  the  satellite,  R  be  the  unit  vector  standing  normally  to  the  orbit 
plane  and 


Q  =  R  X  P  . 

Let  P' ,  0'  ,  R'  be  the  corresponding  vectors  of  the  orbit  of  the  Moon. 
Putting 


Ai(a) 


+1 

0 

0 

0 

+ 

cos 

a 

-  sin 

a 

0 

+ 

sin 

a 

+  cos 

a 

A3(a) 


+ 

cos 

a 

-  sin 

a 

0 

+ 

sin 

a 

+  cos 

a 

0 

0 

0 

+1 

we  have  for  the  components  of  P,  Q,  R  and  P' ,  Q'  ,  R'  in  the  equatorial 
system: 


[P,  Q.  R]  =  AjCn)  •  Ai(i)  •  A3(aj)  , 


(1) 


[P'.  Q'.  R']  =  AjCe)  •A3(n')  •Ai(i')  •A3(aj')  . 


(2) 


The  angle  between  the  equator  and  the  ecliptic  is  e.  Designating  u  as  the 
eccentric  anomaly  of  the  satellite,  we  have  for  the  position  vector 


Let 


P  a  (cos  u  -  e)  +  Q  a  \/l 


e  sin  u 


P  =—  s  +  e'  P' 
a 


(3) 

(4) 


=  P  s  (cos  u  -  e)  +  Q  s  yl  -  sin  u  +  e'  P*  , 

where  s  is  the  parallax. 


s 


a 

I  • 

a 
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Put 


a  ■ 

=  p  •  P'.  =  p  •  O'.  r  =  p  •  R'  . 

The  system  of  Halphen's  formulas  for  the  computation  of  auxiliary  quan 
tities  can  be  slightly  modified  and  rewritten  in  our  notation: 

Kl  = 

_  2  +  e'2  , 

(5) 

K2  = 

(1  -  e'2)(l  -  a2)  -  -  (2  -  e'2)  y2  ^ 

(6) 

II 

y2(l  -  e'2) 

(7) 

g2  = 

i-  (Ki^  -  3K2)  . 

(8) 

K3  = 

(2Ki^  -  9Ki  K2  +  27K3)  > 

(9) 

27  g3 

* 

(10) 

The  next  step  in  Halphen's  method  is  the  computation  of 

°  12'  ’•  1  ■  ^)- 

(11) 

(^) 

-  -ih  H'  2.  1  -f), 

(12) 

and  of 

-yiT  144  / — 

A  =  3  ^  VT  . 

(13) 

J2 

B  =  ^ • 

■^^§2  VS2 

(14) 

The  values  of  f  will  be  proper  fractions.  The  series  (Equations 
12)  possess  an  algebraic  branch  point  at  f  =  0  and  the  convergence  in 

1 1  and 

the 
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The  series  (Equations  15  and  16)  are  alternating  and  the  E-process  can  be 
applied  to  speed  up  their  convergence.  The  formula  of  the  E-process  for 
hypergeometric  series  takes  the  form 


F(a,  b,  c,  -x) 


t  <-■>' 

lt=  0 


(a,  k)  (b,  k)  ^ 
(1,  k)  (c,  k) 


where 


+ 


(a,  N  +  j  ) 
(1,  N  +  j) 


(  b,  N  +  j  ) 
(c,  N  +  j  ) 


N^j  V 

Z_l  2P+1  I  j  }  ' 


(17) 


(q,  k)  =  q(q+l) 


(  q  +  k  -  1) 


Applying  Equation  17  to  the  series  (Equations  15  and  16)  and  putting  N  =  3, 
m  =  19,  we  deduce  the  two  following  expressions  in  which  the  coefficients  are 
rapidly  decreasing  as  the  power  of  w  increases. 


1  ^  f 

-  (  2 

7 

f 

Vl+ 

■  irr 

N/rj 

X  (+  2.3870942 

X  (- 

3.7991784 

-  0.0663082 

W 

+ 

0.3693646 

W 

+  0.0225632 

w2 

- 

0. 1556119 

*2 

-  0.0117691 

+ 

0.0889726 

*3 

+  0.0073743 

- 

0.0586828 

*♦ 

-  0.0051060 

w5 

+ 

0.0419870 

*5 

+  0.0037250 

- 

0.0313364 

*6 

-  0.0027325 

+ 

0.0233758 

*7 

+  0.0019070 

w8 

- 

0.0165247 

*8 

-  0.0011936 

+ 

0.0104483 

*9 

+  0.0006337 

- 

0.0055933 

*10 

-  0.0002710 

+ 

0.0024083 

*11 

+  0.0000884 

- 

0.0007898 

*12 

-  0.0000205 

+ 

0.0001837 

*13 

+  0.0000030 

*!♦ 

- 

0.0000268 

*14 

-  0.0000002 

wi5)  : 

+ 

0.0000018 

*15  ) 
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Tlie  next  step  is  to  compute 


^4  ^  ^3  Kj  K2  , 


and 


K5 

-  Ki  (KjKj  -  3K3)  -  2K2^  . 

“11 

=  K4  {a2-i)  +  K5  +  1  g2K3  , 

^22 

=  1/32  -  1  +  e'2)  +  K  +  A 

62  *^3 

5  2 

1-  e'2 

“33 

=  *^5  ^  g2  [“^(1  -  e 

'2)  +  /32 

^12 

®2  1  =  K4  a/3  , 

®23 

®32  =  (*^4  “  f  82)  /3y . 

«31 

=  «13  =  [^4  “  I  g2  (i  -  e'2)J 

ya  , 

1  ■ 

=  a2  -  1  -  i  K 

^  3  *^1  ' 

®22  • 

:  /32  -  1  +  e'2  -  1  K,  , 

®33  “ 

®12  = 

^21  =  a/3  , 

®23  ~ 

®3  2  “  ^y  ’ 

^31  = 

ajs  =  ya  . 

matrix 

is  formed  with  the  elements 

(1  -  e'2)]  , 


Aij  =  a.  j  A  +  a!  .  B 
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Introducing  the  dyadic  (a  matrix). 


<t 


+  AiiP'P'  +  AijP'O'  +  AijP'R' 


+  AjiO'P'  +  A220'0'  +  A23Q'R' 


+  AjjR'P'  +  A32R'0'  +  AjjR'R' 


(18) 


we  can  represent  the  "disturbing  force,"  averaged  with  respect  to  the  mean 
anomaly  of  the  Moon,  in  the  form 


Fo  =  - 


2km' 


0  •  r 


and  the  averaged  momentum 


=  - 


2km' 

«'3 


p  X  (j)  •  r 


(19) 


(20) 


In  addition,  in  the  system  adopted  in  this  paper,  we  have  to  compute  the 
vector  (Reference  7) 


Ko  =  (l  +f)Fo  -  rr  •  F„  .  (21) 

The  vectors  Mq  and  Kq  must  be  averaged  along  the  orbit  of  the  satellite: 

-Ztt 

=  -^J  Moi-1-.  (22) 


(Koo^ 


I  -2  7r  2  7T 

"  2^  J  *^0  “  2F  I  *^0 

•'ft 


du 


(23) 


This  averaging  is  done  numerically,  by  giving  the  values  0®,  30",  60“  ••*  , 
330“  to  u,  computing  values  of  M^Cr/a)  and  Kplr/a)  for  each  of  these 
angles  and  then  forming  the  arithmetical  means. 


D 
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THE  FORM  OF  VARIATION  OF  CONSTANTS  FOR 
THE  LONG  PERIOD  LUNAR  EFFECTS 


Designating  the  area  "constant"  by  c,  we  have 

c  =  VT  R . 
p  =  a  ( 1  -  e^)  , 


and 


—  =  r  X  F 
dc 


(24) 

(25) 

(26) 


Let  h  be  the  unit  vector  along  the  line  of  nodes  and  let  (i,  j  ,  k)  be  the  unit 
vectors  of  the  equatorial  system.  Then 


h  =  i  cos  fl  +  j  sin  Q, 

R  =  h  X  k  sin  i  +  k  cos  i 
It  follows  from  Equations  Z7  and  28  that 


dR  ,  .  .  dO  „  di 

-  h  sin  1  +  h  X  R  - — 

dt  dt  dt 


and  from  Equation  26  that 


1  •  ■  dn  ,  ,  n  di  e 

h  sin  1  —  +  n  X  R  — — 


dt 


dt 


1-6“^ 


R  4^  +  -^R  ^  =  -i-r 


d  t 


VT 


dt 


X  F 


(27) 

(28) 

(29) 


(30) 


The  long  period  effects  are  deduced  from  the  last  equation  by  applying  the 
double  process  of  averaging  and  taking  into  account  the  fact  that  da/dt  =  0: 


sin  i  =  — —  k  '  Mi 

dt 


'00  ' 


di 

dt 


1 


M, 


00 


VF 

de  _  y/l  - 

d  t  / 

e  a 


h  X  R  , 


R  •  M, 


00 


(31) 

(32) 

(33) 
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It  is  of  interest  to  notice  that  for  long  period  effects,  the  averaged 
momentum  is  the  only  function  of  the  disturbing  force  contained  in  the  Equa¬ 
tions  31  through  33.  Equations  31  and  32  represent  a  particular  case  of 
equations  deduced  by  Makarova  (Reference  8)  for  computing  the  special 
perturbations  of  minor  plcuiets.  From  Reference  7, 

—rR  •  F  -  ^RP  •  K  ,  (34) 

dt  dt  dt  ^  e 

where 


(35) 


We  deduce  from  Equation  34  the  long  period  effect  in  oj,  and  by  taking  Equa¬ 
tion  31  into  consideration  obtain 


P  •  K..  =  ^  j”  P  ■  Ko  ^  ■ 

o 

Equation  36  is  identical  in  form  to  the  equation  for  planetary  perturbations, 
deduced  by  Makarova  on  the  basis  of  different  considerations  from  Musen's 
theory  (Reference  7). 

The  system  (Equations  31,  32,  33  and  36)  is  to  be  integrated  numerically; 
and  the  interval  of  integration  can  be  of  the  order  of  several  days. 


THE  LUNAR  PERTURBATIONS  OF  A  CLOSE  SATELLITE  AND 
THE  SOLAR  PERTURBATIONS  OF  ALL  SATELLITES 

The  computation  of  the  lunar  perturbations  of  a  close  satellite  and  of 
the  solar  perturbations  for  all  satellites  can  be  accomplished  by  developing 
the  disturbing  function  into  series  of  Legendre  polynomials.  For  this  pur¬ 
pose  the  discussion  shall  be  limited  to  the  second  Legendre  polynomial,  and 
the  disturbing  function  takes  the  form 
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„  km'r^/3  9„  l\  /o-t> 

R  =  — jr  cos  H-jj,  (37) 

where  H  is  the  angle  between  r  and  r'  .  Let  u®  be  the  unit  vector  in  the 
direction  of  r'  .  We  have  u®  with  the  accuracy  up  to  e'  in  the  periodic 
terms: 

u®  =  P'  (-e*  +  cos  g'  +  e'  cos  2g')  +  Q'  (sin  g'  +  e'  sin  2g'),  (38) 

and  Equation  37  becomes 

R  "  •  (39) 

The  "disturbing  force"  to  be  used  in  connection  with  the  variation  of 
elements  will  be  ■ 


F 


grad  R 


km'  / a' \ ^ 


(3  u®u® 


r  -  p) 


and  the  momentum  of  F  takes  the  form 


M 


u»u« 


p  . 


(40) 


(41) 


Q 


Substituting 


p  =  P  r  cos  f  +  Qr  sin  f 


into  Equation  41  ,  yields: 

3ktn'  [  a'\^  f—  n  n  n  n 

[P-  u»u».P^  (1  +  cos  2f) 


M 


+  p  X  U®u®  •  0  ^  sin  2f  +  0  X  u"u®  •  P  sin  2f 


+  (J  X  u®u®  .  0  ^  (1  -  cos  2f)j  . 


(42) 


T  aking 


-2tr 


^  f  ^  =  >  *  I 


^  j  ^  cos  2f  dg  =  +  I  e2 
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we  obtain  the  averaged  value  of 

Wo  =  (77)^  [(1  +  4e2)  P  X  u»uO  .  P  +  (1-  e2)  Q  x  u«uO  •  q]  .  (43) 

Let  0,  >Ai  and  6  be  angles  which  the  basic  vectors  P,  0,  and  R  form  with 
the  direction  to  the  Moon;  hence 

cos  0  =  P  •  u®  ,  cos  0  =  Q  •  u®  ,  cos  0  =  R  •  u®  .  (44) 

Replacing  Mqq  in  Equations  31,  32  and  33  by  Equation  43,  and  considering 


h  .  P  X  u® 

h  •  0  X  u® 

(h  X  R)  •  (P  X  u®) 
(h  X  R)  .  (0  X  u®) 
R  •  P  X  u® 
R  -  0  X  u® 
h 


=  +  cos  6  sin  CO  , 

-  +  cos  0  cos  CO  , 

-  +  cos  0  cos  CO  , 

=  “  cos  0  sin  CO  , 

=  +  cos  0  , 

=  -  cos  0  , 

-  P  cos  CO  -  Q  sin  co  , 


we  obtain 

dfl 


sin  1 


d  t 


3km 


2a 


-(7)’ 


di  _  3km' / a' 

~  2a'  ^  \  r'  / 


cos  0 


VF 


de  10  r. - 2  X  m  a"  /a  ^  ,  t  a-j\ 

77= - ^eVl-e  —  — y-pj  cos  0  cos  0  .  (471 


^(l  +  4e2)  COS0  sinoi  +  (l-e^)  cos0  coscoj  ; 


(1  +  4e2)  cos  0  cos  CO  -  (  1  -  )  cos  0  sin  1 


(45) 


(46) 
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'=2  /„.\3 


k  m' / a 


-1041 


-1041 


2C 
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Setting  the  mass  of  the  earth  equal  to  1,  we  have 

k 

-  =  na  , 

and  the  system  (Equations  45  through  47)  becomes 


.  .  dfl 

^  dt 


2v/r 


(l+4e^)  cos  cp  sincti  +  (1-e^)  cosi//  cos  co  cos  &  ;  ^  j 


di 

dt 


-3!!LiL_  [(l  +  4e2)  COS0 


2-Ji 


cosa>—  (1  —  e^)  cosi/fsinct^  cos^; 


n  CO 


(49) 


de  15  m' ne  y 1  -  /  a  /a'\^  ,  , 

dT  - - 2 -  j  j  p  cos 

In  addition  we  have,  with  sufficient  accuracy, 

=  (1-e'^)  ^  +  3e'  cos  g' ; 


and  for  doj/dt  an  equation  analogous  to  Equation  36, 


do)  _  Vp  n  fc-  • 
dT  -  "  —  P  *^0  -  ^  dF 


where 


(50) 


(51). 


(52) 


F  =  ^  (3  u®u®  .  r  -  r) 

r'  ^ 


Substituting 
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15 


2 7e  -  2e^ 

cos‘0)  - — 

2(l-e2) 


It  follows  from  Equations  55  and  56  that 

The  next  step  is  to  compute  the  integral 

i  P-  rr  .  F  d*  =  ^  [  [3(»"  '  «•  I  -Is 

After  some  transformations,  and  taking  Equation  53  into  account,  we  have 

.2tt 

j_  r  ± 

277  J  rp 


(59) 


»  .  rr  •  F  dg 

km'  a 

(l-e2)  r'3 

/  9  2^ 

1  J  cos  p  + 

J  cos^p  -  1^ 

^2‘rr 

.2  77  2 

j  ^  cos 

f  dg  +  (  4-  cos^p 

-  1  cos^p) 

^  j  ^  cos  3f  dg 

Jo  “ 

\  ^ 

^  / 

Jo  ® 

(60) 


From  Equations  57  and  58  and  from 


i-f” 

277  J 


we  arrive  at  a  final  result: 


,277 


COS  3f  dg 


=  - 


^  j  P  •  rr  .  F  dg 

0 

y  C0S^(P  -  Y  cos^i/;  +  2^e  ^  cos^c^  +  y  cos^i/zj  .  ^^^  ) 


km '  i 


( 1  -  )  r 


2^  .<3 


Substituting  Equations  59  and  61  into  Equation  51  and  eliminating  k  in  favor 
of  a  and  n  yields 

dcj  3m'na^  /  a'\^  ,  2j,  2/  /i  2  ■ 

yy  =  — — jy-  ^  j  (4  COS^tp  -  COS^p  -  1)  yj  1  -  B  “  COS!  yy  .  (62) 


The  system  of  Equations  48  through  50  and  Equation  62  should  be  used 
in  investigations  connected  with  the  lunar  perturbations  for  close  satellites 
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and  with  the  solar  perturbations  for  all  satellites.  This  system  contains  all 
the  significant  long  period  terms.  In  addition  to  the  long  period  terms,  the 
lunar  part  also  contains  the  terms  having  periods  of  the  order  of  one  month. 
The  interval  of  integration  for  the  solar  perturbations  can  be  taken  to  be 
approximately  one  month.  However,  for  the  lunar  perturbations  it  would  be 
preferable  to  take  the  interval  to  be  of  the  order  of  days,  to  obtain  a  smoother 
curve  and  to  avoid  ambiguity  in  interpreting  the  results. 


CONCLUSION 


The  methods  described  here  are  more  accurate  and  more  adaptable  to 
the  use  of  electronic  computers  than  the  analytical  development  obtained 
previously.  These  methods  can  be  used  in  investigations  connected  with  the 
stability  of  orbits  and  for  the  separation  of  known  and  unknown  long  period 
effects  in  observations.  It  is  recommended  that  in  practical  use  the  analyti¬ 
cal  development  be  replaced  by  the  semi-analytical  solution  given  here. 


ACKNOWLEDGMENTS 


The  author  takes  this  opportunity  to  express  his  gratitude  ’to  his  col¬ 
leagues,  Mrs.  Ann  Bailie  for  the  careful  preparation  of  the  manuscript,  and 
Messrs.  Arthur  Smith,  Jr,,  and  Lloyd  Carpenter,  who  did  all  the  actual  pro¬ 
gramming  and  the  computations,  thus  making  the  practical  application 
pos  sible. 


-1041 


-1041 


17 


REFERENCES 


1.  Gauss,  K.  F. ,  "Determinatio  Attractionis  quam  in  Punctum  Quodvis 
Positionis  Datae  Exerceret  Planeta  si  Eius  Massa  per  Totam  Orbitam 
Ratione  Tempore  quo  Singulae  Partes  Discribuntur  Uniformiter  Esset 
Dispertita,"  Coll.  Works  3:331-335,  Gottingen:  Konigliche  Gessellschaft 
Wissenschaften,  1876 

2.  Kozai,  Y. ,  "On  the  Effect  of  the  Sun  and  the  Moon  upon  the  Motion  of  a 
Close  Earth  Satellite,"  Smithsonian  Institute  Astrophysical  Observatory 
Spec.  Report  No.  22:7-10,  March  20,  1959 

3.  Musen,  P.,  Bailie,  A.,  and  Upton,  E. ,  "De  velopnnent  of  the  Lunar  and 
Solar  Perturbations  in  the  Motion  of  an  Artificial  Satellite,"  NASA  Tech¬ 
nical  Note  D-494,  January  1961 

4.  Halphen,  G.  H. ,  "Traite  des  Functions  Elliptiques  et  de  leurs  Applica¬ 
tions,"  Vol.  2,  Paris:  Gauthier- Villars  et  Fils,  1888 

5.  Goriachev,  N.  N.,  "On  the  Method  of  Halphen  of  the  Computation  of 
Secular  Perturbations,"  (Russ.)  University  of  Tomsk,  1937 

6.  Goursat,  Edouard,  "Sur  L'Equation  Differ entielle  Lineaire  qui  admet 
pour  Integrate  la  Serie  Hypergeometrique  ,"  Paris:  Ecole  Normal 
Superior  Annal,  2nd  Series  ,  10,  Supplement,  3-142,  1881 

7.  Musen,  P. ,  "Special  Perturbations  of  the  Vectorial  Elements,"  Astronom. 
J.  59(7):262-267,  August  1954 

8.  Makarova,  E.  N.,  "On  the  Simultaneous  Determination  of  Systematic 
Errors  of  Stellar  Catalogues  and  of  Masses  of  Planets  from  the  Ob¬ 
servations  of  Asteroides  ,"  Akad.  Nauk  SSSR,  Biulletin  Instituta 
Teoreticheskoi  Astronomii  7(84):1-18,  1958 

9.  Cayley,  A.,  "Tables  of  the  Developments  of  Functions  in  the  Theory  of 
Elliptic  Motion,"  Mem.  Roy.  Astron.  Soc.  29:191-306,  1861 


MAU  -  Luclej  n*M.  Tk.  0-1041 


-s  o  >5  . 

Jog* 

S  r--  S  ir 

I  is:- 

5  C  O  e 
Z  O  ;3  ft, 

■=s  i3  *•  s 

.2  «  .a  a 

S  <  s  -5 

w  <D  s  ?  i 


c:  H  M  < 
o  U  H 
-  w  3 

2  H 
^  ^ 
Sws 
■S  «  <  I 
6<“ 
•o  iJ  J 

<  o  < 

«'?§ 

o  «  y  ( 

aiS' 

'"Sh 
^  : 
«§<! 

i52?' 

i  s  6>. : 

^  <9  O 

o  c  o  ^ 

'T  ?  Z  Z 

Q  ioO 
H  rt  w  9 : 

^  c  S  2 

o  H  ,1 

W5  «  W  . 
<nZ.X 
Z  Z  O  H  f 


H  E  TO  : 


i  •  -c  o  ' 

s  3  .a  ^  , 

O'^  ** 

:  0/  a>  k  Q>  . 

I  ^  -c  • 


^  *®  o  ! 

•2-“  °  S-i 

3  H  5  *5  j 

i  •  "s  “  I 

”3  111 

,  ^  fci  2  I 


'  t>  -  2:  ft,  I 

"  S  "  S  1 

*Sj  (V  B9  3  < 
t®  C  C  i 

fl»  C 

£  C  "  a  ! 

w  o  ®  §  ' 

o  >v£  ft,  ; 
Q>  .'S  O  JS  I 

b£  —  ^  ' 

*?.  ,D  w  «  I 

£  2 

?  to  a>  g  I 

O  A,  c  u 

C  a  O  a> 

Z  £  r  ad 


B 

I  O  M  ^  f  2T 

2|| 

i  2*  «•«'-  S|« 

rS  ..2  I  8  f -g^ 

‘  "g  Si'S  5  »5^ 

is.sssa|s§ 

;  5  s  O  0  -g  . .  a  «  s 

:s«s«^SS*cav 

I  u  §3^5  vgac* 

!-S3*|S*'C  «w«» 

92.3<0  *^*0  *  * 

“S  <*  C3(Dt« 


"  Ik.  ^  c 

I“i2  ‘ 

e  ^  CO 


lii=s 

iiill 

°s3i»a 

siiSzB 

Z  Z  O  H  A. 


*i- j’! 

shl^l 

u  O  u  ^  « 


JS  Sj._ 

*1  --gjil  « 

£2'»  S^-  S 

3?  s  b"!:!  i: 

iH|5:^ 

3s|«sl- 
g  £  a  2  -2  g  -g 

®  ’S  «♦  *”  o  !!  *• 

jj-g^-'s  £  = 

**  tf  S  k  b 

2  I  =  «:  a  S  £ 

H  E  to  3  £  E  3 


^  2 
«<  ®  *  1  J! 
o2£  It 

S^lii- 

2i aSl  if 

3  ®  t-s  "If 

5 1 "  ^  B  S  ^ 

S  «  a<  ^  ft*  ^9 

i  ;  3  f  I  « 


j  «  I  B  I  '*  I 


2  g  wO  S’ 
•s-  o  -S  -^  2  . 
_§  o  w  «  *5  ^ 

C,  «.  0*  L.' 

.S  t-  5  < 
2  ^  *0  2  -w  • 
"c  ..  c  « 

2  i  §  :  a  I 
2  §  -2 1  s  ^ 

—  i;  <«  p  «  J 

.2  «  a  a  a : 

.tj  <  c  •-  c  d 

c  3  5.  3  « 

w  E  ft*  S  .! 


I  s 

(ft  E 
MO.  ^ 

H  Sc 

(ft  ~  (ft 

z  V  5  ^ 

e  c  . 

o  ft*  « 

Ss  S£ 

.•a«*i5 

OJ  -  W  m 
t«  >  (ft  w 
ft*  (ft  a  (ft 
•—  C  CO  c/3 

A  t*  *0  «  » 

(ft  c  <0  r* 

V  (ft  ^  ^ 


..  c 

1  g  g)0  g> 

»  Q  2  2  ®  ’’ 

u  »  V  u  ' 

5  t-  «J  rj“  .S  i 

2  TI  S  «B  I 

^  ..  c  ^  c  ' 

>.  (ft  —  (ft 

2  i  s :  a  i 

2  §  =  I S  -i 

—  ir  (ft  c  (ft  ' 

£  «  a  a  .a : 

js  <  c  •-  c  i 

p  c  S  c  I 

■S'  CO  E  ft*  E  .! 


5  §¥ 

:  -s  ft* 

d  S.  §  wT 

g  W  a  (;; 

*  t  a  = 

3  (V  7  (ft 

:  Z  4* 

s  e  ^  ' 

9  0>  ,2  e  a* 

*  CM  S  B 

:  «  A.  S 

*  »  .2r  y  al 
2  ki  >  (ft  S 

*  ft*  (ft  a  eft 

5  •-  C  W  CO 
.  b«  ^  « 

jeftCcor- 
j  w  eft  M*  ■<»• 


j  H  .  U 

5<ii 

52  .S 

5  0=® 

5  y  11  o 

jP  Sm 
“1 9  s  ^ 
C  s  *7  < 

"  w  S  3 
5  Z  ft*  5 
3  H  A  5. 


s  3  .a  >  £ 

*-  O'JS  A. 
ft*  ft*  ^  ft*  ^ 

1-2  SH’S' 

£  °  c  -S  . 

ft*  (ft  O  (ft 

5  H  £  i  2; 

I  •  "S  ?  i  d 

S  2  =  a  -g  i 
.2  °  t:  «  “  i 

^  ft*  (ft  o  ®  I 

ft  S  ft*  u  2  1 
w  s  t*  ft*  5;  I 

£  ^  (ft  g  g  ] 
*01  (ft  n  3  c*  ' 
08  p  c  «  ' 

ft,  B 
«  ft*  ft*  ft* 

5  B  -5;  O  U  , 

£  w  >*  a  c  ' 
o  *  5  O  i 

®  >*£  0/  £  1 
ft*  Z  o  .c  S 
«*—  ja  "  ' 
*0  ^  0,  08  , 

£  rt  o  5  .a  * 

^  w  ■>  E  a 

O  Cl® 

5I-2&5 


■5^25-g. 

■“■o  S  2  .£ 

ft*  ft*  .S  : 

5  S  rt  S  g‘ 
o  3|  i-ol 

»  15  a*  §  * 
■o  c  t:  -a  §  » 
£  £  eft  o 
“S  B  c  ”  ®  ’ 

a  2  (ft  T?  ft* 

*  *2  ^  ft*  js 
a  3  o  2  H  • 
c*  u  c  ^ 

0,  ft*  o  a  «  ’ 

£;  a:2  * 

(ft  -o  ^  c  £  ' 

«2  e  35^ 

M  ^  ft*  AA  ft* 

I  a.5 
i  s>st  ” 

®  S  >>  z  . 

*—  o  c  «n  (ft  , 

O  «  I  -S  g 

M  5  s.t3•£■ 

§  "s  a "« s 

Serfs 


e«  ^  " 

I  £  c  £  3 

H  £w5£ 


’  W  Z  2.  ( 
*2d  £  c  j 

'oj  eft  ft9  : 

®  S  « 

Q,  S 

«  a  £  *( 

£  ^  ^  J 

d  ^  ®  ■ 
°  >•£  i 

ft*  z  O  J 

bc—  -O  *■ 

“O  r»  t 
£  £  ^  : 
^  »  00  j 


fi 

a  I  o 
o  c  fefiO 

5  o£  . 

-  |u  .. 

u  S  ^  ^ 

5-1  ■  ^  5 

4»  Q  ^  e  bA 

®*  Z  <  §■  "  c 
.H  So  B  ■ 

S  <  2  i  °  S 
1^3  2s=|i 

S  2  .3  ■  S  S.i 

a  ■<  g  a  e 
•  B  9a: 

M  9  9>  S  ?  £ 


-5  £ 

-  R  §  if 

I  ^  3*  t) 

^  >  S  e 

:Z  S*5^ 
lafo  I  S 


\  a  >  a  Z 

!  a>  cf  a.  n 
i  *2  *3  OT  CO 
•>  >«  *0  ^  « 
I  («  c  (D  r* 

>  U  M  ^  ^ 


lis'Sf  I 

E  o  1 2*15  -ifS 

5  S  5  ^  o  5  ■S 
*  P  9  ttJ  S  —  5 
•—  iJ  <•  g  «  <5  5 


.S  •  £  o 

B  3  .2  f  ° 

tt  ^  Si  5 
S  ^  5  E  H  rt  ' 

r  hi  A  ^  a; 

a>  o  A  5  ^3 
S  ^  a>  A  o  n  I 

!  5  S  .s  S  £  i 

f  i  ?  i  i 

-  “I  Is  11 

y  «  O  h  "S 

H  «  £!  «  I  i!  : 

«  9  ^  B  U  J 

-  ^  A  »  3  t)  S 

5  «  C  C  A  * 

0  4>  .  9  c 

=  IC”sg' 

5  “  o  «  g  H  j 

i  °  o  £  - 

a>  iS  o  .c  SJ  t 
5  bc—  ,Q  —  *  a 
i  ?  .B  w  CO  *  *Q 

i  5  A  O  5  .2 

•  >  aS  «  £ 

g  0,  g  5  g  a 

*  5  .c  .2  a>  2  « 

>  ^  5  iT  a.  £  s 


|S0  ;ii|S35 

S-S-sl^S® 


C3  H  U  < 
o  O  H 
9  W 
A  bi  **"* 

is  s  ■ 

.2  w  S 

.S(5<| 
6<” 
•a  J  j  ' 
<  O  < 

U  M  ^  C 

s-iS' 

T3  a ' 
S  Q  3  S 
o  „  s 

S25  ! 
-.igfe  j 

A*  A  ” 

O  C  o  ^ 

^  O  y  Z 
-I  6  r?  c 
Q  a>  O  2  a 

z<  j  P  : 

H  «  w  9  5 

.  c  ac  s  . 

^  o  a  ^ 

CO  ;s  ^  W  i 

<  A  Z  PC  a 

Z  Z  O  H  a 


O  S 

|i? 

ItL 

e  «>>  • 

A-  o  4)  « 

Sr  65 
*”*  Si  ^  9 

h  >  A  5 
«  A  a  A 
*2  B  CO  CO 
U  *0  » 

A  B  (O  t- 
U  A  A<  A* 


2  S  -a  .  £ 

O  U  2  »-  S  J!  S.- 
isSS  °-2s  =  53 

2i*<Ort  I 

“‘-"“S.Si"! 

=  .i“|j 

--f*  stiaii  “* 
a  £  S  3  g  ,2  -o  I 
0  **  o  3  Ob'S  *®i 

2g;S 

=  1  ®  35  £  g  2  S 

«<Z-AtS-M 
I  .  —  ^  A  a  3  O  C 

*»^0  B9CCA*' 

s  5  *  S  Q>  §•-  a»  2 

I  9££'*^2!,a>  2*^ 
'^’S««6£ 

I  ’3  T3  «fc-  O  3  "  C 

!  I  g  °  £■£  <U  £  * 

!  ^  »  S>55£  "  a 

:  SS'S-g'siS  £  * 
"3|“gEg3 
E  SiSsa  g.5  § 


S 

s  t  o 

O  e  b£0 

i  I  .2  -■ 

Jog® 
'Z  «•  a» 

•w  t-  «  jr- 
.2  ID  X 

•o  ..  C  ® 
^  >.  A  r 

^iij 

-iii 

.3  3  u  I, 
9  <  c  -3 
c  9  2 
9-  (o'  £  «  i 


2  9  •- 
-  o  9 

:  -2  <u 

8  M  a  ^ 

3  e 

3  A  7  S 

Iz  ^  2^ 

t  c  . 
j  oT  o  2  ? 
^  csi  r  s  5 


ti<  >  A  r 
«  A  a  A 

-  C  w  CO 
.  t«  "O  ... 
A  C  (O  t^ 
y  A  A*  A* 


s 

,o  g  ^<5  ^ 

J  (5  ’§  2*  “ 


J  (5  S  00 15 

l"i§i 

^  ^  A  r  A 

S9  I  §  c 

5  S  2  9  2 

^  p  r  a>  r 

•2  «  .2  §.  2 
.ta  <  9  •-  c 

B  9  9,  3 

ir  (cf  s  y  £ , 


>  5  c 

5  -2  g 

**-  a  I  "»* 

9  A  ^  2*  y 
2  9  'S 

2  "o  5  ?■  « 

g‘  I  s'  S  II 

if  £ “II 

:::  w  V  «  a  CO 
«  Z  •-  e  CO  oo 
to  .  u  'O  .  . 
C  CO  A  C  CO  t^ 
.—  CA  y  A  A*  A* 


£  ■“  J,  y  . 
y  .o  5  ■« 

4,1  ”  S: 
^  ZS  a 

*-  3  A  £ 

1^  A  (>  l- 
o  y  jS  y  • 
a> 

•o  c  t;  *0 
S  .3  «  o 

9  f-  ••^  ' 


£  2  o  «i 
“■  l.  c 

£  a5-2 

"  T3  ^  C 

«  .2  E  3  i 

—  (U  _  • 

=  S-c  E  • 

P  V 

S  c -S  « 

O  9 

5  9  «  ( 
®  2  ®  w  ! 

«  5  S  i  «. 

S  'a  y  *“  1 
y  o  ^  aj 


o  2  . 
>  -2  9  y 

H  E55 


:  •«  o  A  ®  ^ 
^  -5  E  3  .3  »  5 

!  -9  u  o*  JS 

e-3".gEH« 

cj  t  2  ^ 

9  5  o  -S  5  'A 

052^»«  C’O 

ill 

.  -9  _  U  A  •"*  9 

2  ■;  4; !;  3  4; 

—  3'5't;—  firS 
"-m*2“e3: 

»'r*7c'a‘*®p4, 
*”  *<  o  Abba 
o  A  y  _  9 

^  A  y  y  "a  y 

n  p-sis « °  g. 

£|^ro«5E 

w  £  9  O  y 

«  ?  "  4,  si  S  "  ■ 

5cwbCsB^** 
o  y  ^  .  .  MS  w  ' 

^  ^  —  5a^5.2. 

^^*5a  «"9  9' 

2o-£gyg5g 

2  6  5«5r  S.5 


-I  w 

f-  . 

d  y 

u  H  y 

o  w  g. 
a  £  Q. 

s<  2 

S”S=. 

5  <  A* 

0  M  n’  ® 

5  E  Q 

3  Z  "  ^ 
>  < 

^  fc,  3  < 

'o  *’9 
■  z  •  ^ 

)  y  y  O 

ip  5w 

jis^ 

•  s  .  < 


p  £  , 

>.1;  B  , 
y  .0  5  ■« 

•0*8. 
y  y  B 
■9  A  —  '2 

9  A  E 

b«  A  u  ki 

®  "  r  5  ’ 

A  Tj  "y 
*9  C  t-  ’C 

|a  «  o 


9  9  O  A  I 
9  ^ 

£ 

A  .«  O  C 

2-2SS: 

9  2-9  £  ■ 
g  9>  y 

5  ft  S*'! 
*—  5  9  A 

°  j|  ”, 

E^  4,  =  1 
y  ®  B  y  ! 

A  9  ^  9  5 


g  2  9  y 

H  s5S, 


j  a  g  g  "S  f  o 

i  g  ®l  " "I 

=  r.'a’S-sEH|« 

l.2^  ft-  *2 

iS"£-°°c’3& 

‘yA^yAOAA 

!  p^af  g  =  4, 1 

Ar<*-*c^  —  Afc^ 
BA?  •O  O 

lapS^ja^ESc; 

yaSSlfi  CE’D  c'a 
.■S.«^‘"P""S22 
•  y®  o  fcA**5-2 


B«i<»'yAc9’2 

i-niooigii 

fe2  *SOAAc 

I  O  A  (U  9  B  ‘ON 

2  B  O  ^B  y  9  5 

®  A  5*r  y  t.  ■"  C 
>"ysaoB3yy 
;cAM“.D'^**yB 
g5‘gS.-AAy  9 
r  —  5  ®i3.2.-  o 

S-sl-gSo^'S 

ia|2§£®SS 

BnZrrars'C 


